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ABSTRACT 

We report the discovery of a linear relationship between the root-mean-square (rms) 
variability amplitude and the mean flux in the accreting white dwarf binary system 
MV Lyrae. Our lightcurve, obtained with the Kepler satellite, spans 633 days with 
quasi-continuous 58.8 second cadence resolution. We show, for the first time, how 
this cataclysmic variable displays linear rms-flux relations similar to those observed in 
many other black hole binaries, neutron star binaries and Active Galactic Nuclei. The 
phenomenological similarity between the rms-flux relation observed here and in other 
X-ray binaries suggests a common physical origin for the broad-band variability, inde- 
pendent of source type, mass or size of the compact accretor. Furthermore, we infer the 
viscosity parameter, a, and disk scale height, H/R, using two independent methods. 
In both cases, both values are found to be uncomfortably high to be accommodated 
by the disk instability model. 

Key words: accretion discs, stars: binaries: close, stars: novae, cataclysmic variables, 
stars: MV Lyrae, black hole physics, stars: oscillations. 
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Compact interacting binaries (CBs) are close binary systems 
usually consisting of a late-type star transferring material 
■ onto a compact object such as a black hole (BH), neutron 
star (NS) or a white dwarf (WD) via Roche-lobe overflow. 
With an orbital period on the order of hours, the trans- 
ferred material from the donor star forms an accretion disk 
surrounding the compact object in order to conserve angular 
momentum. The dynamics and physics governing the flow 
of matter accreting onto compact objects is however still 
debated. BH and NS compact binary systems are powerful 
and highly variable sources in X-rays, and it is by studying 
the fluctuations in X-ray luminosities that we have made 
progress toward s understand i ng the physics invo lved within 
accretion disks ( Belloni 2 0071 : Ivan der Klis!l200"7T ) . For exam- 
ple, thanks to the Rossi X-Ray Timing Explorer (RXTE), 
Ivan der Klis et all (| 1996T ) have been able to detect fast (kHz) 
quasi-periodic oscillations (QPOs) in NS binaries, associat- 
ing these with the inner edges of the accretion disk. The 
detected kHz QPOs can imme diately be used to constrain 
the mass and radius of the NS (|Miller et alj|l998h . 
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Variability is a common characteristic of X-ray binary 
systems (XRBs, consisting of NS or BH accretor), cata- 
clysmic variables (CVs, consisting of a WD accretor) and 
Active Galactic Nuclei (AGN, where the accretor is a super- 
massive black hole). Specifically, the power spectrum of this 
variability consists of two components: single or multiple 
quasi-periodic oscillations(QPOs), and an aperiodic broad- 
band noise component spanning several decades in temporal 
frequency. QPOs have been extensively studied in a num- 
ber o f accreting sources, ra nging from WD binaries (jWarnerl 
|2004 iPretorius et alJ|2006h. to super-massi ve black holes at 
the centre of AGN |Gierlinski et al.l l2008h . On the other, 
hand the properties of the broad-band noise component has 
been examined in detail only in XRBs and AGN. 

One feature of the broad-band noise X-ray variabil- 
ity which seems to be ubiquitous in XRBs and AGN is 
the so-called rms-flux relation, which is a linear relation- 
ship between the absolute root-mean-square (rms) ampli- 
tude of variability and the flux, suc h that the sources become 
more variabl e as th e y get brighter dUttlev fc McHardvll200ll; 
Uttlev et all 120051 : Icieissner et al l 120041 : iHeil fc Vaughanl 



2010h The rms-flux relation is remarkably linear and ap 



plies over a broad range of time-scales, so that the rms of 
high frequency variations correlates with those at low fre- 
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quencies l|Uttlev et alj|2005h . Crucially, the rms-flux rela- 
tion is observed while the power-spectrum shape is itself 
stationary, implying that the rms-flux relation is a funda- 
mental property of the variability process, and not connected 
with other processes which can change the variability prop- 
erties on longer time-scales. Finally, it can be shown that 
the rms-flux relation predicts a log-normal distribution of 
fluxes and implies that the variability process is non-linear 
l)Uttlev et al.l l2005h . The current most-favoured model to 
explain the rms-flux relation, as well as many other fea- 
tures of the variability i n XRBs and AGN, is the propagat- 
ing fluctuations model (|Lvubarskii|[l997l : iKotov et aLll200ll : 
lArevalo fe Uttlevll2006h . In this model, the observed vari- 
ability is associated with modulations in the effective vis- 
cosity of the accretion flow at different radii. More specifi- 
cally, the model assumes that, at each radius, the viscosity 
- and hence accretion rate - fluctuates on a local viscous 
time scale around the mean accretion rate, whose value is 
set by what is passed inwards (again on the viscous time 
scale) from larger radii. The overall variability of the accre- 
tion rate in the innermost disk regions is therefore effectively 
the product of all the fluctuations produced at all radii. 

For CVs, the presence of broad-band noise in the power 
spectrum has been recognised for many years and is usu- 
ally referred to as "flickering" in the literature. However, 
the methods employed to study XRB lightcurves have not 
been applied to CV data. The properties of flickering have 
been mainly studied throug h the use of eclipse maps (eg. 
iBaptista fc Bortolettol l|2004h L which have identified the ex- 
istence of two different and independent sources of flickering: 
the low and high-frequency flickering components. The low- 
frequency component originates from the overflowing gas 
stream from the secondary star, and is associated with the 
mass transfer process. The high-frequency component, on 
the other hand, originates from the inner accretion disk re- 
gion, and no evidence of variable emission from the hot spot, 
gas stream, or white dwarf ha s been observed (c ontrary to 
some previous suggestions, see lBruchlll99ll2000h . 

XRBs and AGN, which do show the rms-flux relation, 
share many properties with CVs. Both XRBs and CVs are 
observed to exhibit large-amplitude outbursts, which are 
theoretically explained as arising due to thermal/ viscous in- 
stabilitie s in the accretion disks s urrounding the compact 
objects (|Shakura fc Sunvaevlll973l ). Furthermore, similarly 
to XRBs, CVs also display outflowing jets in conjunction 
with specific change s in the spectral energy distribution 
l|Kording et al . 2008). Therefore, it is natural to suspect that 
CVs should also show behaviour associated with a fluctuat- 
ing accretion flow, including the rms-flux relation which is 
seen in accreting BH and NS. However, although CVs share 
many similarities with XRBs and AGN, they also differ sub- 
stantially in other respects. Specifically, the radius of a WD 
is over an order of magnitude larger than that of a NS or 
the event horizon of a BH. Because of this, accreting white 
dwarfs are not subject to the effects of strong gravity. Addi- 
tionally, accreting WDs have not been observed to display a 
comptonised component as in XRBs and AGN, and further- 
more they do not posses the extreme magnetic fields which 
are found in some NS. Therefore it remains an important 
open question as to whether accreting WDs should show 
the same rms-flux behaviour as seen in XRBs and AGN. 

The broad-band noise component, present in XRBs and 



AGN, is well studied because of the high throughput X- 
ray missions with good timing capability (RXTE, XMM- 
Newton), allowing very short cadence (down to ~ms) obser- 
vations, at high sensitivities. In the optical. iGandhi (120091) 
obser ved three XRBs with ULTRACAM jPhillon etafl 
I2007D and, as expected, also found that the lightcurves dis- 
play similar rms-flux relations as those observed in the X- 
rays (although the optical emission could possibly originate 
from the jet and not the disk). Accreting white dwarfs are 
faint in X-rays, since their gravitational potential wells are 
not as deep as those of XRBs (so they show lower inner 
disc temperatures) and also (unlike AGN and XRBs) they 
do not seem to possess strong Comptonising coronae. Be- 
cause of this, detailed timing studies in the X-ray region 
have not been possible for disk-fed white dwarfs, and most 
timing studies of these systems concentrated on optical light. 
Specifically, a series of papers (Warner 2004-; IPretorius et al.l 
2006) already noted analogies between the QPOs observed 
in CVs and those observed in XRBs. However, detailed anal- 
ysis of the broad-band noise component in CVs was so far 
not as straight forward as for XRBs. The reason for this 
is twofold: first, the broad-band noise component displays, 
in the optical, low-amplitude, few per cent fractional rms 
variations requiring sensitive cameras; second, the broad- 
band noise component can only be studied with long, con- 
tinuous, observations. This is because optical broad-band 
noise is thought to be generated further out in the accretion 
disk compared to X-rays, and is thus observed on longer 
timescales. Furthermore, broad-band noise is a stochastic 
process, requiring statistical studies of many realisations. 
This, in turn, demands long continuous observing. In order 
to truly assess whether accreting white dwarfs display the 
rms-flux relation, and thus determine if this phenomenon 
is indeed ubiquitous in all accreting compact objects, irre- 
spective of accretor type, a sensitive, long, high-cadence and 
continuous light curve in the optical light is required. This 
kind of observation ha s only now becom e possible, thanks 
to the Kepler satellite dKoch et alj|2010h . 



Here we analyse the Kepler light curve of the accreting 
white dwarf system MV Lyrae (one of about 14 CVs known 
to exist within the Kepler field-of-view) , and show that it 
does indeed exhibit a n rms-flux relation sim ilar to that ob- 
served in Cyg X-l by iGleissner et all (|2004h . MV Lyrae is 
a VY Scl nova-like system, spending most of its time in the 
high state (V f» 12 — 13) and occasionally (every few years) 
undergoing short-d uration (weeks/mon ths) drops in bright- 
ness (V w 16 - 18, iHoard et al]|2004h . It is not clear why 
MV Lyrae, and all VY Scl type stars, undergo this kind 
of behaviour, but it is known that these systems have ex- 
tremely low ma ss transfer rates at their minimum brig htness 
(6 x 1O" 1O M iLinnell et ai1l2005l ; IHoard et all 12004 ). The 
inferred distance to MV Lyrae is m 530pc, with a secondary 
star cooler than 3500K (|Hoard et alj|2004l ). The orbital pe- 
riod of this CV is 3 . 19 ho urs, with a binary mass ratio of 
0.43 (|Skillman et alj|l995l ). We observed MV Lyrae during 
its final ascent, in a low-to- high luminosity transition, where 
its optical lig ht mostly originates f rom its nearly face-on ac- 
cretion disk l|Skillman et al.| [l995). In the next section wc 
will describe the data acquisition and analysis, whilst our 
results and discussion are presented in Section [3] 
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2 DATA ANALYSIS 



The MV Lyrae light curve was obtained with the Kepler 
satellite and provided to us by the Science Operations Cen- 
trcQ in reduced and calibrated form aft er being run throug h 
the standard data reduction pipeline (I Jenkins et al.l I2010T) . 
The response function for the photometry covers the wave- 
length range 4000-9000A. Here, we only consider the Single 
Aperture Photometry (SAP) light curve and do not attempt 
to analyse the Pre-search Data Conditioning (PDC) light 
curve, as the latter is optimised for planet transit searches. 
As a consequence of this choice, data artifacts which would 
have been removed by the PDC, remain in our Kepler SAP 
light curve. These artifacts are the result of rapid attitude 
tweaks to correct pointing drift, caused by the loss of fine 
pointing accuracy. This pointing accuracy loss is due to the 
spacecraft thrusters firing regularly every 3 days to desat- 
urate the reaction wheels and can affect Kepler data for 
several minutes. Furthermore, data gaps occasionally occur 
due to Kepler entering anomalous safe modes. During such 
events no data is recorded and following data is often cor- 
related with electronic warming. Further details of artifacts 
within Kepler light curves can be found in the Kepler Data 
Release Notes 3jj. Here, we make no attempt to correct these 
artifacts, but simply remove them from the light curve. 

FigQj, shows t he short cadence (58.8 seconds; 
iGilliland et all 1201(f ). barycentre corrected, light curve 
for MV Lyrae obtained during the first four quarters of 
Kepler operations. The lightcurve spans an interval of 633 
days. The visible data gaps present in the light curve are 
due to the artifacts described above as well as the monthly 
data down-links. At first glance the light curve reveals an 
obvious steady increase in flux (over 1 magnitude) and also 
an increase in the flux variability (rms) during the rise. 

To characterise the relationship between the rms vari- 
ability and the mean flux, we follow ed a similar proce- 
Gleissner etafl (|2004l ) on Cyg X-l. Whilst 



dure as used bvinieissner et all [|^UU4| | on Uyg A-i. wnust 
IGleissner et alj (2004) used the mean PSD level to measure 



rms, we directly measure it from the spread in the lightcurve 
datapoints. First, the MV Lyrae lightcurve was split into 10 
day segments. For each segment we computed 10 minute 
flux averages and corresponding rms. All 10 day segments 
analysed displayed linear rms-flux relations after binning 
the 10 minute averages by flux. One example is shown in 
Fig[2] We carried out the same exercise employing smaller 
and larger data segments, and longer and shorter averaging 
timescales. In each case, we still recovered linear rms-flux 
relations. This is already strong evidence that the optical 
variability observed in the accreting WD is analogous to the 
X-ray variability seen in XRBs and AGN. 

In order to quantify the linear relati onship in Fig. wc 
mode lled it with a function of the form (jUttlev fc McHardvl 
[2003), 

a = k(F-C) (1) 

where a is the rms value, F is the flux in e~s _1 , whilst 
k and C are the fit parameters. This linear model fits the 
general shape of the majority of our flux-sorted segments 

1 The data is propritary. The first five quarters of Kepler obser- 
vations will become public on January 7 2012 

2 http://archive.stsci.edu/kepler/documents.html 
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Figure 2. Rms-flux relation obtained from a 10 day segment 
(indicated by the yellow shade in FigfTJ in MV Lyrae. 



reasonably well, although occasionally large x 2 values are 
recovered, implying there exists some intrinsic scatter. This 
scatter can be caused by any non-stationarity in the variabil- 
ity processes which causes variations in the rms-flux relation 
on the 10 day time-scale of the segments. Furthermore, the 
presence of any frequency-varying QPO could break the lin- 
ear rms-flux relation by varying the power spectrum over the 
frequency-range used to measure the rms. In Fig(3]we show 
the noise-subtracted power spectral density (PSD) of the 
four coloured, independent, sections shown in the lightcurve 
of Fig[T] The timescales used to measure the rms are those 
at frequencies higher than 10 minutes, marked the dashed 
black line in Fig[3] For each segment, these rms values will 
respond to flux variations on longer time-scales, down to the 
10 day length of the segment. 

The PSDs shown in Fig(3] show possible candidate 
QPOs at low frequencies and br eaks just above 10~ 3 Hz . 
From the modelling performed bv lArevalo fc Uttlevl £2006), 
these breaks can be associated with the viscous timescale at 
the inner-most edges of the accretion disk. Above the break 
frequency, no annuli contribute significant variability power, 
and the PSD bends downwards. Given that the four different 
quadrants appear to display shifts in the break frequency, it 
is possible that mass accretion rate variations are displac- 
ing the inner edges of the accretion disk. Furthermore, at 
frequencies below the break, there appears to be candidate 
QPOs. These features at low frequencies, and the PSD break 
shifts, will be analysed and discussed in a paper in prepa- 
ration: here we are only concerned with the high-frequency 
component observed in MV Lyrae. 

Another interesting feature observed in MV Lyrae con- 
cerns the evolution of the rms-flux relation with both flux 
and time. In Fig(4]we show the best-fit lines for the rms-flux 
relations in all of our 10-day segments. In order to give a 
feel for the temporal evolution of the rms-flux relation, each 
colour represents a different section in the lightcurve, with 
the colours in Fig[T]associated to those in Fig(3] Specifically, 
there appears to be a hysteresis-like behaviour in the evo- 
lution of the rms-flux relation within the broad-band noise 
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Figure 1. MV Lyrae light curve obtained by the Kepler satellite. The source was observed in short cadence mode (58.8 second cadence) 
for 633 consecutive days. The light curve segment within the yellow shaded region was used to produce the rms vs. flux plot in FigfS] 
The blue, red, green and magenta regions have been used to produce the PSDs and rms-flux relations in Fig[3j Fig|4]and Fig[5j The last 
233 days, marked by the magenta region, has been used to produce the flux distribution in Fig|6] 




Frequency (Hz) 



Figure 3. Power spectral distribution obtained from the MV 
Lyrae lighcurve shown in Fig[T] The colour coding is analogous 
to that in Figfl] so that each independent PSD was obtained by 
the corresponding light curve segment. The Poisson noise level, 
ranging between 2 — 3 X 10 -5 , has been subtracted. In this paper 
we analyse frequencies higher than those marked by the black 
dashed line. 

of MV Lyrae. After about 150 days, the rms-flux relations 
seem to change behaviour and detach from the original trend 
(marked with blue lines in Fig(4]) and move to higher fluxes 
maintaining similar rms variability. In proximity of the peak 



flux in Fig[T] after about 250 days of observations, the rms- 
flux relations seem to change behaviour once again during 
their flux descent (red lines), seemingly looping back to their 
original track after having reached a maximum rms variabil- 
ity of about 6%. This behaviour has not been observed in 
any other accreting compact binary or AGN, and it is not 
clear at the time of writing what the physics behind this 
phenomenon is. In Fig[5]we show the gradient and intercept 
for the best-fit lines in Fig[4] colour coded in the same way. 
Si milarly to the re s ults o btained on the BH binary Cyg X-l 
bv lGleissner et all |2004:) , the gradient and intercept in MV 
Lyrae seem to be correlated as well. This suggests that some 
fundamental and analogous process is giving rise to the rms- 
flux behaviour in both XRBs and CVs. Interesting is also the 
fact that in the gradient vs. intercept plane there seem to be 
two distinct "tracks": the rms-flux relations obtained from 
the blue and red lightcurve segments seem to have a sys- 
tematic offset in the gradient when compared to the green 
and magenta segments. Although we cannot properly quan- 
tify this behaviour at the moment with the available data, 
the evolution in Fig0] also seems to suggest that the rms- 
flux behaviour systematically changes during the lightcurve 
evolution. 

Finally, we also examined the flux distribution within 
the MV Lyrae lightcurve. Specifically, the observed fluxes 
in XRBs and AGN which do show the rms-flux relation 
are log-normally distributed. In order to perform this test, 
we produced a histogram of fluxes from the last 233 days 
(marked by the magenta region in Fig[l} of the MV Lyrae 
lightcurve. We then fitted the histogram with both a sym- 
metric Gaussian and a three-parameter log-normal model. 
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Figure 4. Right: The best-fit lines for the rms-flux relations in all of our 10-day segments, where one is shown in Fig. [2] In order to 
follow the temporal evolution of the rms-flux relations, colours are matched to those in Fig[T] The black data points are the mean flux 
and mean rms values of the computed rms-flux relation fits. Left: Same figure, where the y-axis displays percentage fractional rms. 
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Figure 5. Best-fit parameters (gradient and intercept from Eq[T] 
normalised by the mean flux of each 10-day segment) from the 
rms-flux relations shown in Fig|4] The colour coding is matched 
to those in Fig[T] 
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Figure 6. Binned frequency distribution of the mean source flux 
for the magenta section of the MV Lyrae light curve in Fig[T] 
(black data points). A log- normal distribution fit is shown in red 
and a Gaussian fit in blue. A log-normal fit statistically describes 
the data better than a Gaussian. 



The number of flux measurements per bin (N) was required 
to be at least 50 for the fits, with bin error assigned as 
vN. We performed this exercise on the last 233 days of the 
lightcurve, as these appear to be most the stationary in mean 
flux: a requ irement for the log-normal behav iour to be ex- 
pected (see lUttlev et alj 12003 : iGandhil l200sh . Statistically, 
the Gaussian fit resulted in the goodness-of-fit parameter 
X 2 /d.o.f of 18802/269, whilst the log-normal fit resulted in 
334/269. MV Lyrae clearly exhibits a log-normal flux distri- 
bution. The excess residual of our fit can be explained by the 
weak non-stationar ity o f the ligh t curve , as also explained in 
lUttlev et all (|2005l ) and lGandhil (|2009l ). 



3 DISCUSSION 

We showed how the accreting WD MV Lyrae obeys a lin- 
ear rms-flux relation like many other BH/NS compact in- 
teracting binaries and AGN. This is the first time such a 
correlation was found in an accreting object other than a 
BH/NS X-ray binary or AGN. Furthermore, we showed how 
the rms-flux relation evolves with time, displaying, what ap- 
pears to be, a hysteresis-like behaviour. The physics of this 
phenomenon is not yet understood, and it is the first time 
it has been observed in an accreting compact object. Addi- 
tionally, we confirmed our findings by ascertaining that the 
flux distribution of MV Lyrae follows a log-normal distribu- 
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tion, as seen in other accreting compact objects displaying 
the linear rms-flux relation. 

The leading interpretation of the rms-flux rel ation is 
that of a fluctuating accretion disk (|Lvubarskrilll997l ). where 
variations in the accretion rate occur at all radii and then 
propagate inwards on the viscous timescale within the ac- 
cretion disk. These fluctuations couple together in a multi- 
plicative way as they propagate down to the centre, so that 
low frequency fluctuations produced at large radii modulate 
the higher frequency fluctuations produced further in. The 
multiplicative nature of the process will yield flux variations 
which follow a log-normal distribution, like the one in Fig[S] 

IChurazov et al.l (|200ll ) have pointed out that, in a stan- 
dard thin disk, viscous damping will smooth accretion rate 
fluc tuations on the v iscous timescale. This means that, for 
the lLvubarskiH [|l997h model to work, the fluctuations must 
happen in a geometrically thick accretion disk flow (such 
as an advection dominated flow or thick disk). We can use 
the PSD break just above 10 -3 Hz (ss 13 minutes, shown in 
FigGJ, or the relative rms variations, to constrain the vis- 
cosity parameter, a, and the disk scale height, H/R. 

The PSD break in FigO can be associated to the vis- 
cous timescale at the inner edges of the accretion disk 
(|Arevalo fc Uttlevl l2006). If this were the case, we would re- 
quire a geometrically thick disk of about H/R ~ 0.1 (where 
H/R is the ratio of the vertical scale height to radius in 
the disk), and a viscosity parameter a w 1. This is shown 
in Fig(7] where we computed the viscou s timescales result- 
ing from the standard a-disk solution of iFrank et al.l (|2002T ) 
(Eq. 5.68) as a function of H/R for an accretion disk ex- 
tending all the way down to the WD surface at 7 x 10 8 cm. 
Both these values are considered to be high for a standard 
dis k solution, and are pushing into uncomfortable limits of 
the lShakura fc Sunvaevl (|l973l ) parameter space. 

Following the p rocedure adopted by 

IBaptista fc Bortolettol (120041 ) on V2051 Oph, we also 
used the rms variability amplitude to obtain a and H/R. 
From the inferred 3% rms variability (see Fig(4| we obtain, 
similarly to V2051 Oph, a « 0.16. This estimate assumes 
an H/R value of about 10~ 2 , and is based on the effects of 
magneto-hydrodynamic (MHD) turbulence in an a c cretio n 
disk, as modelled by ICeertsema fc Achterberel (l992). 
However, when the variability reaches 6% in MV Lyrae, 
a ss 0.62, again pu s hing i nto uncomfortable limits of the 
IShakura fc Sunvaevl (|l973T ) parameter space. To make a 
consi stent with the expected viscosity of hot disks (|Lasotal 
|20Qll), the disk scale height, H/R, would have to increase 
to 0.1. This would then be in line with a geometrically 
thick, hot, disk. If this were the case, the rms variability 
amplitude would be an indicator of the disk scale height 
(assuming a to remain constant at it's hot value of w 0.1 
- 0.2). More importantly, all the above estimates for the 
disk scale height and a are suggesting that the propagating 
fluctuation models may have difficulty in explaining the 
variability if the disc behaves like a "standard" thin 
disc. One way to overcome these inconsistencies would 
be to invoke an overflowing gas stream, as mentioned by 
IBaptista fc Bortolettol (|2004h . This flow could exist above 
the standard thin accretion disk, extending all the way 
close to the WD surface, and possibly provide the thick 
disk required to explain the fast viscous timescale inferred 
from the PSD breaks. 
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Figure 7. Viscous ti mescale as a functio n of H/R from the stan- 
dard a-disk solution (Frank et al. 2002j, Eq. 5.68). Each colour 
corresponds to a different a with black, blue, green and red cor- 
responding to 0.001, 0.01, 0.1 and 1 respectively. The horizontal 
line marks the PSD break at 13 minutes. 



One other important implication of our discovery re- 
gards the interpretation of flickering noise in CVs in general. 
It has been recognised for many decades that almost all CVs 
display at least some optical flickering. Many suggestions 
have emerged in order to explain the origin of this flickering, 
by associating it to the accreti on hot spot or the boundary 
layer c lose to the WD surface |BrucrJll992l . l200d ). More re- 
cently, IBaptista fc Bortolettol (|2004h showed that CV flick- 
ering is caused by the existence of two different and inde- 
pendent sources. The low-frequency flickering component is 
caused by the overflowing gas stream from the secondary 
star, whilst the high-frequency component by the inner-most 
regions of the accretion disk. Our result is in agreement with 
the interpretation that the high-frequency flickering is cre- 
ated within the accretion disk. The fact that we observe the 
linear rms-flux relation within MV Lyrae (a canonical flick- 
ering CV) gives us confidence that the fluctuating accretion 
disk model is also at play within other CV systems. Fur- 
thermore, the observed fractional rms amplitude of ~ 3%, 
generally observed in MV Ly rae (see FigBj), is also similar t o 
that found in V2051 Oph by IBapti sta & BortolettcJ l |2004h . 

Although it is clear from the iLvubarskiil (| 19971 ) model 
that the variability is the result of modulations from an en- 
semble of coupled perturbations, the physical mechanism of 
variable emission is not yet known. Some suggested that syn- 
chrotron emissi on at the base o f the outflowing jet could be 
a possibility l|King et al.l 120041 ). However MV Lyrae is not 
known to display a jet (or synchrotron emission), and it is 
not likely that a comptonised component is driving the rms- 
flux relation out to the large radii observed in the optical. 
More likely instead is that instabilities within the accretion 
disk are responsible for driving the continuum variability. 

Thi s interpretation is a lso supported by the recent re- 
sults of lUttlev et all l201ll . They showed that for XRBs 
the blackbody component in X-rays, associated with the 
disk, substantially leads the power law component associ- 
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ated with the comptonised region close to the compact ob- 
ject. Thus, the comptonised region, thought to drive jet- 
launching, is responding to changes in the accretion disk, 
and not vice-versa. In this respect, we would expect MV 
Lyrae to display energy-d ependent time lags in optical light 
l|Arevalo fc Uttlevl l|2006l )). where the broad-band variability 
arising from bluer bands is lagging the variability observed 
from the redder ones: simply a consequence of the fact that 
the fluctuations are propagating inwards within the disk, 
affecting the redder outer edges before the bluer inner ones. 

The rms-flux relation remains an enigmatic observa- 
tional feature of all accreting compact objects, but it clearly 
contains information about the dynamics of the infalling ma- 
terial. Although the microphysics of accretion disk dynam- 
ics are not yet fully understood, our result strongly suggests 
that the processes involved are the same in all accreting 
compact objects, independent of mass, type and size. Our 
result has thus revealed yet another connection between ac- 
creting compact objects, suggesting that the physics of disk 
accretion are universal. 
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